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Abstract 

A molecular tool that can compare the abundances of different DNA sequences is necessary for compar- 
ing intergenic or interspecific gene expression. We devised and verified such a tool using a quantitative 
competitive polymerase chain reaction approach. For this approach, we adapted a competitor array, an ar- 
tificially made plasmid DNA in which all the competitor templates for the target DNAs are arranged with a 
defined ratio, and melting analysis for allele quantitation for accurate quantitation of the fractional ratios 
of competitively amplified DNAs. Assays on two sets of DNA mixtures with explicitly known compositional 
structures of the test sequences were performed. The resultant average relative errors of 0.059 and 0.02 1 
emphasize the highly accurate nature of this method. Furthermore, the method's capability of obtaining 
biological data is demonstrated by the fact that it can illustrate the tissue-specific quantitative expression 
signatures of the three housekeeping genes G6pdx, Ubc, and Rps2 7 by using the forms of the relative abun- 
dances of their transcripts, and the differential preferences of \gf2 enhancers for each of the multiple \gf2 
promoters for the transcription. 
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1 . Introduction 

Gene transcription levels have been measured by 
methods based on Northern analysis or reverse tran- 
scription (RT)- polymerase chain reaction (PCR) 
assays. These approaches are generally effective for 
comparing the RNA levels of any gene among many 
samples. However, these methods are unsuitable for 
com paring the quantitation data generated indifferent 
places or times, primarily due to the paucity of quanti- 
tation standardization. Comparing the expression 
levels of different genes is even more problematic, 
even though this task is important for understanding 
the coordinative nature of gene expressions and the 
relative strength of regulatory elements.' Microarray 



technology and recent advancements in sequencing 
technology, such as RNAseq, have opened the way for 
a thorough examination of the expression of many 
RNAs; they also facilitate the visualization of relative ex- 
pression levels of multiple genes.^'^ Nonetheless, these 
assays are yet to be tailored to attain sufficient accuracy 
or reproducibility between different assay platforms; 
they are also unsuitable for examining a small 
number of genes in numerous samples.'*'^ 

The quantitative competitive PCR (qcPCR) method 
measures the amount of target DNA of interest relative 
to the amount of an artificially made DNA competitor. 
With this measurement, the absolute quantity of 
the target DNA can be determined if the amount 
of the corresponding competitor is known. ^'^ 
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A single-nucleotide difference between two otherwise 
identical DMAs is the best competition structure that 
warrants the equal amplification efficiency of these di- 
morphic DMAs in the PCR.The resultant PCRamplicons 
have been analysed for the allele ratios in the samples 
using a wide range of application platforms including 
Maldi-TOF, pyrosequencing, and real-time DNA 
melting.^"^ ' However, qcPCR approach needs to deal 
with some inevitable uncertainties. The quantitation 
of the competitor templates or samples of interest is 
largely performed in a spectrophotometric manner; 
that is, the method takes into account the common 
measurement biases that occur between measure- 
ments and between laboratories. 

In this study, we present an alternative method of 
nucleic acid quantitation based on a qcPCR approach. 
With this method, the amount of a given sequence in 
a mixture such as cDNA or genomic DNA (gDNA) is 
determined in relation to the amount of another se- 
quence in the mixture. To accomplish this task, we 
introduced a competitor array as a novel reaction 
component in the qcPCR. The competitor array is an 
artificially made plasmid DNA in which the competi- 
tor templates for the target DNAs in question are 
arranged with a defined aspect ratio. This structure 
of the competitor array enables one to infer the quan- 
titative relation between two or more DNA sequences 
if the relative quantity of each sequence to its respect- 
ive competitor templates in the competitor array is 
determined. We also exploited high accuracy and pre- 
cision of the melting analysis for allele quantitation 
(MAAQ) method in determining the relative quantity 
of the target DNA with respect to the competitor tem- 
plate. The MAAQ method utilizes the differential 
melting patterns of either the fluorescent resonance 
energy transfer (FRET) probe encompassing on a 
single-nucleotide variation or restriction fragment 
length polymorphic DNAs. This method mathematic- 
ally analyses the melting curve of a sample to extract 
the constitutional fraction with which each allelic 
version contributes to the curve. It offers a 
simple and rapid procedure as it does not require 
post-PCR purification or additional biochemical reac- 
tions while ensuring high accuracy and precision. In 
addition, the MAAQ analysis can be performed with 
a relatively less expensive device than other 
methods. Various exemplary assays are presented to 
emphasize the validity and versatility of this method 
in the field of nucleic acid quantitation. 

2. Materials and methods 

2.1. PI asm ids 

The plasmid DNAs containing a single target or 
competitor template were made by TA cloning of 
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the corresponding PGR amplicons into pCR2.1 
(Invitrogen) so that the inserted fragments could be 
liberated with EcoR\ restriction. Competitor template 
was generated by PCR amplification by applying a 
mismatch-carrying primer. Plasmid DNAs containing 
multiple template sequences which were used for an 
equimolar mixture of the target and competitor tem- 
plate and for the competitor arrays were constructed 
by a means of a single or repeated subcloning 
process. For the subcloning, a Spe\/Xba\ fragment of 
a plasmid was placed into the Xba\ or Spe\ site of 
another plasmid. Plasmids that contained any form 
of an inverted repeat were avoided. Further informa- 
tion on the plasmids used in this study is presented 
in Supplementary Table S5. 

2.2. RT and PCR 

RT was carried out on 1 |jLg of the total RNA with an 
iScript™ cDNA synthesis kit (BioRad) in accordance 
with the manufacturer's instructions. PCR was per- 
formed with ~0.2 jjlI of the reverse transcribed cDNA 
template, 20 ng of gDNA or 0.1 pg of plasmid DNA 
for appropriate analyses and 0.5 jjlM of corresponding 
primer pairs: G6pd-1 F and G6pd-7R for G6pdx; Ubc- 
1 F and Ubc-4R for Ubc; Rps-4F and Rps-2R for Rps27; 
Igf2-E1F1 and Igf2-E4R1 for the Igf2 PI transcript; 
Igf2-E2F1 and Igf2-E4R1 for the Igf2 P2 transcript; 
Igf2-E3F1 and Igf2-E4R1 for the Igf2 P3 transcript; 
and Igf2-E4F1 (GTACCAATGGGGATCCCAGT) and Igf2- 
E4R1 for the Igf2 total transcript. Information on the 
primer sequences is presented in Supplementary 
Table S5. The competitor arrays were prepared by 
cleaving the corresponding plasmid and subsequently 
diluting it in water containing Xphage DNA (1 ng/ixl). 
For the competitive PCR, each template (except for 
the references) was mixed with its appropriate com- 
petitor array and amplified by PCR by using 40 cycles 
at 94°C for 1 0 s, 56°C for 2 0 s, and 72°C for 3 0 s. 

2.3. Melting analysis 

The relative ratios of the target DNA to its competitor 
in the PCR amplicon were determined as previously 
described. Briefly, 5 |jlI of PCR amplicons were mixed 
with 5 |xl of 1 0 mM ethylenediaminetetraacetic acid 
and with 0.4 |xM of FRET probes, Igf2-Sl and Igf2-Al 
(see Fig. 1 b) for Igf2; G6pdx-S, Red640-5'-GCACTG 
TTGGTGGAAGATGTCA-3', and G6pdx-A, 5'-CACACG 
GATGACCAGCTCATTACG-3'-FITC, for G6pdx; Ubc-S, 
Red 64 0- 5'-TACGTCTGTCTTCCTGTTTGAC- 3', and 
Ubc-A, 5'-GGCTCTTTTTAGATACTGTGGTGAGG-3'- 
FITC, for Ubc; and Rps2 7-A, 5'-CCTTCTGATGA 
ATGTGGTGCTGGAG-3'-FITC and Rps2 7-S, Red640-5'- 
TCTGGGAAGCCACTTTGACAG-3', for Rps27. The DNA 
melting was performed with Roche LightCycler 2.0 
under the following conditions: 30 s at 95°C; and 
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Arrangement of the compclilor templates for A, 
B and C sequences in a plasmid with a given 
aspect ratio of a;P:y to form a competitor array. 



Liberation of each competitor template in the 
competitor array by endonuclease reaction 



Addition of the processed competitor array into 
a DNA mixture such as cDNA that contains the 
target DNA molecules with an unknown ratio 



J Co-amplification of the competitor template and 
its target for each sequence 



Determining the relative quantity of a target to 
its compclilor via the MAAQ method 



Determining the relative ratio of the target 
sequences in the mixture by reflecting the aspect 
ratio of competitor template in the competitor 
array 




Expected 



Figure 1. Determination of relative abundances of different DNA sequences, (a). A schematic representation of the entire procedure for 
determining the relative abundances among three sequences (A, B, and C) In a cDNA sample. The respective competitor templates 
(A', B', and C) which are arranged In a competitor array with a given aspect ratio of a:fi:y are liberated by endonuclease restriction, 
mixed with cDNAs, and subjected to PCR reactions for each sequence of Interest. The relative quantities of the target sequences with 
respect to the corresponding competitor In each ampllcon, n, p, and q, as assessed by IVIAAQ, are fed into a calculation based on the 
aspect ratio of the competitor template in the competitor array to yield the relative abundances among the target sequences, (b). A 
set of artificial DNA mixtures constituted by the restriction of the plasmids (A-J) that contain cDNA fragments of three 
housekeeping genes with variable copy numbers as denoted under the name of the mixtures and fractionated on an agarose gel. 
The identity of each DNA Is denoted on the left. (c). The relative abundances of two cDNA sequences were determined and scatter- 
plotted against the expected values on the log scale axes. 



8 min at either 45°C (for /g/S.cept/x, and L'fac)or44°C 
(for Rps27) at 0.1 °C stepwise increments in tempera- 
ture until the probe's melting is completed. The 
background-subtracted fluorescence values (640/ 
Back- 5 30) from the Red 640 dye were extracted with 
the aid of LightCycler Software 4.05 and then 



submitted for calculation as previously reported.^ ^ 
The relative abundance of two sequences in each 
sample was determined from the ratio of the T/C 
values of the two comparing sequences after reflecting 
the copy numbers of their corresponding competitor 
templates in the competitor array used. 
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3. Results 

3.1 . Determining the relative abundance of DNAs in 
simple mixtures 

Figure la depicts the typical procedure for deter- 
mining the relative abundances of the sequences of 
interest in a mixture. To demonstrate the feasibility 
of this quantitation approach, we used it to determine 
the relative quantity of a cDNA sequence with respect 
to another sequence in the artificially made cDNA 
mixtures. A set of cDNA mixtures was composed 
through the construction of 10 plasmids (A-J); they 
contained variable numbers of each of the three 
cDNA fragments of mouse Ubc, Rps27, and G6pdx 
genes. Plasmids H and I were from independent 
clones with the same sequence composition. The rela- 
tive abundances of the sequences in each plasmid are 
therefore explicitly known (Fig. 1 b). Another plasmid 
DNA, the competitor array CAl , was made to 
contain a single copy of each competitor template, 
namely a single-nucleotide substituted version of the 
respective cDNA sequence. These plasmid DNAs 
were structured to liberate each cDNA fragment as a 
single entity upon EcoR\ digestion (Fig. 1 b). Each 
mixture generated by the EcoR\ digestion of the corre- 
sponding plasmid was individually mixed with the 
£coRI-restricted CAl and subsequently subjected to 
competitive PCR reactions along with three standard 
templates, namely the target cDNA, the correspond- 
ing competitor template, and an equimolar mixture 
of the target cDNA and its competitor template. 
One of the fluorescently labelled FRET probes for 
each cDNA sequence was designed to span the nu- 
cleotide of variation between two competing tem- 
plates so that the probe's melting from the 
complementary sequences are differentiated depend- 
ing on the completeness of the base pairing. The 
melting curves that we generated from the PCR ampli- 
cons by applying the respective FRET probes were fed 
into the MAAQ process to determine the relative 
quantities of each target cDNA sequence with 
respect to the competitor template (T/C for short). 

Supplementary Table SI presents the T/C values of 
each cDNA sequence (G6pdx/C, Ubc/C, and Rps2 7/ 
C) in the cDNA mixtures; the quantities were deter- 
mined from the quadruplicate PCR reactions. The 
abundances of a cDNA sequence relative to another 
{Rps27 /GGpdx, G6pdx/Ubc, and Ubc/Rps27) in each 
mixture were then obtained from the ratios of the re- 
spective T/C values of the two sequences in question. 
Figure 1 c shows a scattered plot of the observed 
versus expected values of the relative abundances; 
the plot highlights the accuracy of this assay. The 
average relative error for this assay is 0.059, demon- 
strating that our qcPCR approach can accurately 
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assess the abundance of a sequence in relation to 
another sequence in a simple mixture. 

3.2. Determining tlie relative abundance of different 
DNAs in complex mixtures 

We subsequently considered whether our approach 
would be reliable in the analysis of complex samples 
such as those that are likely to be encountered in a 
laboratory. To address this matter, we measured 
the copy number ratios of the autosomal Ubc and 
the X-linked CGpdx sequences in male and female 
gDNAs. We chose gDNA because it is a particularly 
complex mixture of sequences and provides a strin- 
gent test for our assay. The Ubc/G6pdx ratios were 
expected to be 1 and 2 for females and males, re- 
spectively. About 20 ng of each liver gDNA of five 
females and five males was mixed with 0.05 pg of 
CAl and subjected to PCR reactions for two 
sequences. We then applied the MAAQ. Five XY 
samples have an average ratio of 1.978, and five XX 
samples have a ratio of 1 .000 with an average relative 
error of 0.02 1 . These results strongly confirm that this 
method can measure the relative abundances of the 
target DNAs, even in a highly complex mixture such 
as gDNA (Table 1). 

3.3. Determining the tissue-specific expression 
signatures via the forms of relative abundance of 
cDNAs 

Because of the proven accuracy and reliability of our 
method, we used it to obtain two types of biological 
data. The first type of data is the relative expression 
abundances among three housekeeping genes in dif- 
ferent mouse tissues. A preliminary assessment using 
the serial dilutions of the CAl on the tissue cDNAs 
indicates that the abundances of the Ubc and Rps27 
transcripts are in an equivalent range whereas the 
abundance of GGpdx transcript is far less than those 
of the two transcripts (more than a hundredfold in 
some cases, data not shown). To accurately cover 
this wide range of expression levels, we made an add- 
itional competitor array, CA2, to carry the copies of 
GGpdx, Ubc, and Rps2 /sequences at a ratio of 1:7:7. 
The PCR amplicons of each sequence were generated 
by using a mixture of each cDNA and CA2 as a tem- 
plate. Subsequent analysis revealed the T/C values 
and the relative abundances among the target 
sequences (Supplementary Table S2). Figure 2a visua- 
lizes the relative abundances of two subjected 
sequences in the cDNAs of the liver, heart, and brain 
of three female adult mice. These results demonstrate 
the tissue-specific expression signatures of these three 
genes with forms of their relative transcript abun- 
dances; they also demonstrate the relatively small 
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Table 1 . The relative abundance of nucleic acid species can be reliably calculated in gDNA samples 



Female Male 





F1 


F2 


F3 


F4 


F5 


Ml 


M2 


M3 


M4 


M5 


G6pdx 






















7G6pdx 


0.601 


0.51 2 


0.51 8 


0.510 


0.543 


0.378 


0.385 


0.21 4 


0.326 


0.400 


SD 


0.008 


0.007 


0.010 


0.01 3 


0.01 0 


0.020 


0.01 5 


0.01 1 


0.01 2 


0.01 5 


G6pdx/CA1 


1.504 


1.048 


1.074 


1.039 


1 .1 87 


0.607 


0.625 


0.273 


0.484 


0.666 


Ubc 






















7Ubc 


0.599 


0.51 1 


0.519 


0.506 


0.548 


0.561 


0.553 


0.336 


0.49 


0.567 


SD 


0.008 


0.007 


0.01 7 


0.008 


0.006 


0.006 


0.010 


0.023 


0.01 1 


0.014 


Ubc/CAl 


1.492 


1 .044 


1.081 


1.024 


1 .21 3 


1.276 


1.238 


0.506 


0.960 


1 .31 0 


Ubc/cepdx^" 


0.992 


0.996 


1.006 


0.985 


1 .022 


2.1 02 


1.981 


1 .857 


1.985 


1.957 


Obs/Exp 


0.992 


0.996 


1.006 


0.985 


1 .022 


1.051 


0.990 


0.929 


0.992 


0.984 



A sample of gDNA was prepared from the kidneys of 1 0 adult mice (five females; F1 -F5 and five males: Ml -M5). The rela- 
tive abundance of the Ubc and C6pdx genes was then determined as described in the text. The terms -yCSpdx and 7Ubc are 
the proportions of the target sequences within the competitively amplified amplicons, and G6pdx/CA1 and Ubc/CAl are 
the relative ratios of the CGpdx and the Ubc sequences with respect to their corresponding competitors, respectively; the 
standard deviation (SD) is based on from three measurements; the term Obs/Exp refers to the ratios of the observed per 
expected values. 
^Average relative error, 0.021. 
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Figure 2. Relative abundances among the transcripts of three housekeeping genes that display tissue-specific expression signatures, (a) A 
plot of the relative expression ratios of two of the three housekeeping genes in three adult tissue samples as determined from three 
individuals. The error bar denotes the standard deviations among the individuals, (b) A plot of the relative expression ratio as 
determined from multiple RT reactions of a single individual; the plot highlights the precision of the assay. The error bar denotes the 
standard deviations among multiple RT reactions. 



Individual variations [which have a 0.125 coefficient 
of variation (CV)]. 

Figure 2b shows the relative abundances among the 
transcripts of the three tissues of a single Individual. 
The values were averaged from four measurements. 
We determined three of the measurements from 
three Independent RT reactions by adapting another 
version of competitor array, CA3, with the respective 
copy ratio of 1:9:9 for the sequences of G6pdx, Ubc, 
and Rps27. The CV of these measurements (0.079) 
confirms that our quantitation method has a high 
level of precision. 



3.4. Measuring the relative promoter strengtli of 
multiple promoters of a gene 
We tried to obtain the relative transcriptional activ- 
ity of the multiple Igf2 promoters In embryonic tissue. 
The utilization of multiple promoters for transcription 
Is not rare In a eukaryotic system.'^ The usage rate of 
each promoter could be the subject of developmental 
or tissue-specific regulations. If the promoter strength 
Is linearly correlated with the cellular level of the tran- 
script It drives, our assay scheme is well suited for 
assessing the relative strength among the promoters. 
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This type of assessment provides insight into the rela- 
tive promoter activity of the given samples and con- 
firms the potential of our approach as a versatile 
quantitation method. Transcriptional regulation of 
the murine Igf2 gene is relatively well characterized 
with respect to the role and position of c/s- regulatory 
elements, including promoters and enhancers.^ ^ Fetal 
and neonatal Igf2 transcription is derived from mul- 
tiple promoters, PI, P2, and P3, through appropriate 
interactions between regulatory sequences; these 
types of interactions provide a good model system 
for comparing the transcriptions of multiple promo- 
ters and specific interactions.^"^ 

We conducted an assay in which each promoter-spe- 
cific transcription was compared with the total level of 
Igf2 transcription in two embryonic tissues, the liver 
and muscle with endodermal and mesodermal 
origins, respectively. The primers were designed to 
span two consecutive exons to prevent any potential 
noise from gDNA contamination. The promoter speci- 
ficity of the PCR amplicon was given by placing a part of 
the forward primers on the promoter-specific exons 
(Fig. 3a). The competitor templates for the promoter- 
specific amplification were made by placing only a 
single-nucleotide substitution on exon E4; they were 
then arranged into a plasmid to constitute the com- 
petitor array CA4. Each promoter-specific competitor 
was also used as the competitor for the entire Igf2 tran- 
script by applying a new forward PCR primer placed on 
the 5'-end of the shared exon E4 (Fig. 3a). In this way, 
CA4 can be exploited in two different competitions 
with an aspect ratio of 1 :3. One competition is for the 
promoter-specific transcript; the other is for the total 
transcript. Each cDNAwas added to the CA4 and subse- 
quently subjected to four separate PCR amplifications: 
three of the PCR amplifications were specific for the 
promoter-specific transcript and the other was used 
for the total transcript. After determining the T/C 
values of each PCR reaction, we obtained the relative 
abundance of each promoter-specific transcript with 
respect to the total Igf2 transcript; the results are pre- 
sented in Fig. 3 b and Supplementary Table S3. 

The Igf2 transcription in the liver necessitates the 
involvement of an endoderm enhancer, whereas the 
muscle expression largely relies on a mesoderm en- 
hancer. Our results confirm that minimal transcrip- 
tion activity occurs on promoter PI, regardless of 
which enhancer is involved; specifically, only 6-7% 
of the total transcript is attributed to promoter PI ac- 
tivity in both the liver and muscle. On the other hand, 
the two types of tissue with developmentally different 
origins have different preferences with regard to the 
utilization of promoters P2 and P3: compared with 
P3, P2 has more transcription in the muscle but less 
in the liver. This difference indicates that the relative 
strength of the two promoters depends on the type 



of enhancers used. Moreover, the results confirm 
that this method is highly reproducible; the standard 
deviations of the measurements of three independent 
RT reactions are minute (ranging from 0.051 to 
0.001, CV= 0.033). 

If exon E4 is spliced solely for the transcripts derived 
from the three promoters addressed here, the relative 
transcription levels of the three promoters should add 
up to 1. The fact that the aggregate transcription 
levels for the muscle samples are 1.010 and 1.008 
is consistent with this notion. On the other hand, 
the aggregate transcription levels for the liver 
samples are 0.920 and 0.938. This difference may 
be attributed to the possible existence of one or 
more transcripts that contain exon E4 in the liver. 



4. Discussion 

In this study, we presented a qcPCR method that 
can be used to measure the relative abundances 
among the sequences of interest in any DNA 
sample. Because of the technical difficulty of obtain- 
ing a direct quantitative comparison between two 
unlinked sequences, we exploited the absolute relativ- 
ity of their respective competitor templates which are 
arranged in the competitor array. In addition, the high 
accuracy of the MAAQ method was exploited for the 
quantitation of the fractional ratios of the competi- 
tively amplified DMAs. The relative abundances of 
two different sequences were consequently obtained 
in a series of plasmid constructs and gDNAs. These 
abundances turned out to be very accurate; the 
average relative errors of 0.059 and 0.021 , respective- 
ly, confirm the reliability and practical feasibility of 
this approach. 

The quantitation scheme was further validated by 
its successful acquisition of biological data. We com- 
pared three housekeeping genes in terms of the abun- 
dances of their transcripts for multiple types of mouse 
tissue and analysed Igf2 expression in terms of the 
proportion of the promoter-specific transcript in the 
total Igf2 transcripts. These results efficiently illustrate 
the tissue-specific expression signatures of the three 
housekeeping genes and the relative transcriptional 
strength of the Igf2 promoters. 

With an exploitation of the high level of accuracy 
and precision, our method can potentially be used 
as a form of standardized quantitation. That is, the ex- 
pression level of any gene of interest can be deter- 
mined relative to the expression of a single gene or 
a combination of several reference genes in the 
same sample. Thus, our method can produce a 
stand-alone value for the transcription level of a 
gene in a particular sample. Furthermore, this value 
can readily be compared with other results that are 
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Figure 3. Promoter-specific transcriptions as measured in relation to the total transcription level, (a) The fetal lgf2 expression as derived 
from the promoters PI , P2, and P3. Each exon (El , E2, and E3) constitutes a promoter-specific part of tgf2 mRNA that is spliced to the 
shared exon E4. The primers (arrows with dotted line in the middle) were designed to span two consecutive exons except for one 
forward primer, Igf2-E4F, which is placed on the 5'-most part of exon E4 to account for the total transcription of lgf2 and the 
common reverse primer. The detailed structure of the exon E4-exon E5 junction region regarding the positions of the reverse 
primers Igf2-E4R1 and Igf2-E4R2 and the FRET probes Igf2-A and Igf2-S is depicted in the lower part of the figure, (b) The relative 
quantities of the promoter-specific transcription with respect to the total transcription yield the relative abundance of each 
transcript in the samples. The relative transcript levels, which were measured from muscle tissue and liver tissue, represent the 
mesodermal and endodermal origins, respectively. 



independently generated by the same method at any 
place or time. In this regard, Ubc and Rps27 can be 
considered as reference genes because their relative 
transcript abundances were relatively constant in all 
the tested tissues. Our data also demonstrate that 
the assay remains robust across a wide range of ex- 
pression levels (e.g. the G6pdx/Ubc value in liver is 



0.005). This range can be further extended by apply- 
ing a competitor array that contains the respective 
competitor templates with a larger aspect ratio. 

For the purpose of sharing and exchanging the 
quantitation data of the transcripts, various applica- 
tions have been developed to present the quantities 
on an absolute scale such as copy numbers per 
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cell.^^'^^ Miura et al.^^ reported a qcPCR approach 
with which the transcription level of thousands of 
yeast genes were described as copy numbers per 
cell. Kanno et o/. also presented a method to 
convert the microarray data or quantitative PCR data 
to copy numbers per cell. Both approaches utilize a 
series of standard RNAs which is spiked into the 
sample of interest in order to provide a way of calibra- 
tion of inter-experiment differences or normalization 
of quantity data in the experimental sets. Although 
these approaches are effective in providing the struc- 
tural architecture of transcriptome abundances on 
an absolute scale, it still requires standardization of 
standard RNAs in order to cross-refer those data to 
others that are independently generated. In this cir- 
cumstance, it would be more practical to present 
the transcription level of a gene in relation to that of 
another in the same sample as we did. Many 
approaches including above two examples and 
others that are focused on defining copy number vari- 
ation in the genomes are indeed capable of delineat- 
ing the abundances of multiple sequences in the view 
of their relativity.^ ^'^ ^ However, the accuracy and the 
precision that these approaches can confer are not as 
high as those our method can, or were not addressed 
probably due to disregard of presenting the relative 
abundances of intergenic sequences. 

On the other hand, application of our method for 
the standardized quantitation is yet to address 
certain problems that are potentially caused by the 
differences in the cDNA conversion efficiency of indi- 
vidual RNAs and between the segments of a given 
RNA. The different reaction conditions or types of 
enzymes used for RT consequently produce variations 
in the relative quantity of given sequences in the 
cDNA pool. Nonetheless, our approach provides a 
way of examining the extent to which the relative 
abundance of those sequences varies on account of 
the different reaction conditions; it also provides a 
way of determining the best conditions for minimiz- 
ing such variations. 

With simplicity, speed, and accuracy, our quantita- 
tion approach provides a straightforward paradigm 
for comparing gene expression levels across cell 
types, developmental stages, and even different 
species. Most importantly, with this approach, experi- 
ments performed in different laboratories or at differ- 
ent times can be correctly referenced to one another. 

Supplementary Data: Supplementary Data are 
available at www.dnaresearch.oxfordjournals.org. 
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